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We analyse numerically a pinch-type instability in a semi-infinite planar layer of inviscid 
conducting liquid bounded by solid walls and carrying a uniform electric current. Our 
model is as simple as possible but still captures the salient features of the instabil¬ 
ity which otherwise may be obscured by the technical details of more comprehensive 
numerical models and laboratory experiments. Firstly, we show the instability in liquid 
metals, which are relatively poor conductors, differs significantly from the astrophysically- 
relevant Tayler instability. In liquid metals, the instability develops on the magnetic 
response time scale, which depends on the conductivity and is much longer than the 
Alfven time scale, on which the Tayler instability develops in well conducting fluids. 
Secondly, we show that this instability is an edge effect caused by the curvature of 
the magnetic field, and its growth rate is determined by the linear current density and 
independent of the system size. Our results suggest that this instability may affect future 
liquid metal batteries when their size reaches a few meters. 


1. Introduction 

The electromagnetic force that results from the interaction of electric current with its 
own magnetic fie ld is us ually rota tion al and, thus, induces a flow when the conducting 
medium is fluid ( Boiarevics et ail 1989h . However, some rotationally or translationally 
invariant current distributions can produce a purely potential electromagnetic force which 
can be balanced by the pressure gradient alone. Such quiescent equilibrium states are 
not always stable and can collapse when subject to an arbitrary small disturbance. The 
prominent example is the electro magne tic pinch, which can disrupt or kink conductors 
carrying strong electric currents (Haines et al 2000 ). It usually affects conducting media 
with deformable boundaries and has a growth rate determined by the balance of inertia 
and the characteristic magnetic pressure. This results in the instability that de velop s with 
the Alfven wave speed regardless of the electrical conductivity of the fluid ( Tavlen I960l l. 


Electromagnetic pinch-t ype insta bility can develop also in the incompressible liquids 
bounded by solid walls ( Michael!Il954 1. In a perfectl y co nducting liquid, this internal 
pinch instability also develo ps on the Alfven time scale ( Velikhovlll959h . The previous two 
studies were extended by Tavler (1961) to general non-axisymmetric instability modes in 
a per fect ly co nduct ing fluid bounded by solid cylindrical walls. Later it was suggested by 
Vandakurovl ( 197211 and iTavler ( 1973!) that a similar instability can affect the interiors 
of stars containing toroidal magnetic fields. The strong stratification in stellar interiors 
makes this instability nearly horizontal and^hus, significantly different from the pinch 
i nstab ili ties with deformable boundaries ( Tavler 1957 *1. In order to stress this difference, 
Spruit ( 1999 *1 termed the astroph ysic al variety of pinch instability the Tayler instability. 
This term was later used by iRiidiger et al\ ( 2007t f in a much broader sense to refer 
to current-driven instabilities in homogeneous fluids including liquid metals. The latter 





















































2 


J. Priede 


type of instability, which differs from the Tayler instability not only by the absence of 
radial stratification but also by its resistive nature, which will be dis cuss ed in t his p aper, 
was presumably observed in the liquid-metal experiment bv ISeilmaver et~cd. d2f)121. The 
notion of Tayler instability is further broadened in the recent study by Herreman et all 
( 2015 1 who refer by it to all pinch-type instabilitie s st udied by R. J. Tayler including 
also the classical case with deformable boundaries ( Tayler 19571 ll96Clh . It is important 
to note that surface deformation provides an additional pinch instability mechanism. 
Although surface deformation always involves fluid flow, it is not required for the internal 
pinch-type instability, which can be driven by the fluid flow alone. In ideally conducting 
fluid, where the magnetic field is frozen in, flow causes the same-order disturbance of 
the magnetic field as_that asso ciated with the surface deformation, i.e., proportional to 
the displacement ( Tavlerl 1957 1. In a poorly conducting liquid, the perturbation of the 
magnetic field caused by the fluid flow is proportional not to the displacement but to 
the product of velocity and conductivity. This perturbation is much weaker than that 
associated with surface deformation ( Tavlen 1960 1. 

Another peculiarity of internal current-driven instability is its reliance on the curvature 
of the magnetic field. For example, in the liquid metal annulus carrying an axial current 
and bounded by solid walls, instability vanishes when the radius of t he g ap tends to 
infinity and, thus, the circular magnetic flux lines straighten out ( Priedel 20151) . This 
is consistent with the absence of purely magnetic instability in the planar perfectly 
conducting liquid layer permeated by a straigh t co-planar magn etic field whatever its 
distribution over the thickness of the layer ( Qgilvie fc Pringlg 19961) . In contrast to 
cylindrical geometries, where the curvature of the magnetic field is inherent and pinch- 
type instabilities usually occur, in planar geometries, such instabilities are likely to be 
induced by the edges and corners, around which the magnetic field bends. Moreover, 
rectan gular configuration s may be relevant to the recentl y developed liquid metal bat 
teries ( Wang et aLll2014l) . for which the electrode mixing ( Kellev &; Sadow avl 20141 ) due 


to potential current-driven instabilities is one o f the major concerns (IStefani et a/.ll2011 


Weber et al. 20131 2015 : Herreman et al. 20151) . 


In this paper we consider an elementary model of pinch-type instability consisting of 
a semi-infinite planar layer of inviscid conducting liquid that is bounded by solid walls 
and carries a uniform electric current. This is the simplest possible model of the internal 
pinch instability in planar geometry. It allows us to elucidate the basic characteristics 
of the instability, which may otherwise be obscured by technical details of more realistic 
numerical models and laboratory experiments. In particular, we show that the internal 
pinch instability is an edge effect caused by the curvature of the magnetic field, and its 
growth rate is determined by the linear current density and independent of the system 
size. This suggests that the instability can be prevented either by reducing the field 
curvature at the edges or by limiting the layer thickness but not its lateral size when 
the areal current density is fixed. We also point out that the instability is significantly 
different in well conducting and highly resistive fluids like liquid metals. In the latter, 
the instability development time depends on the electrical conductivity rather than being 
determined by the Alfven wave speed as for the Tayler instability. This somewhat limits 
the astrophysical significance of liquid-metal laboratory experiments. On the other hand, 
it implies that this type of instability, if any, in the liquid metal batteries significantly 
differs from the Tayler instability in astrophysics. 

The paper is organized as follows. Mathematical model is introduced in fJ5] and the 
linear stability problem is formulated in fj3] Numerical method and the main results are 
presented in JI] and (}5| The paper is concluded with a summary and brief discussion of 
the results in SJS] 

























































2. Formulation of the problem 

Consider a channel of width L and depth S between two solid vertical parallel walls filled 
with an inviscid liquid of density p and electrical conductivity a. The liquid carries electric 
current of a uniform density j 0 directed upwards along the rigid edge of the channel, 
which coincides with the y-axis of the Cartesian system of co-ordinates shown in figure [T| 
The layer is assumed to be thin with S <C L and, thus, treated as a planar undeformable 
sheet with the areal density p = pS and the respective electrical conductivity a = ad. The 
electric current with the linear density Jq = e y joS that flows along the sheet generates 
a magnetic field whose normal ( z ) component through the sheet can be found using the 
Biot-Savart law as 


B 0 (x) 


fJ-oJo f L dt; 
2tt Jo £-x 


Mo A) 
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In 



^qJq 

2n 



( 2 . 1 ) 


where the last approximate expression holds in the vicinity of edge for x <C L, which 
will be the main focus of this study. On the other hand, the sheet approximation implies 
x S. Note that expression m represents a superposition of the fields of straight wires 
which may be thought to form the sheet. As only the gradient of the field B' 0 (x) « — 
is required in the following, the width L drops out of the model. It means that the layer 
is treated as semi-infinite. The interaction of the electric current with its own magnetic 
field gives rise to the pinch force with the areal density F o = J o x B o = e x JoBo(x), 
which is balanced by the x-component of the effective pressure gradient 


Pq(x) = -JqB 0 (x). 


( 2 . 2 ) 


Note that for sheet, the effective pressure P is defined as force per unit length. 

Disturbance of this quiescent equilibrium state by an arbitrary slow liquid flow, which 
is confined to the plane of the sheet by the solid bounding walls, is governed by the 
linearised 2D Euler’s equation 


pd t v i = -VPi -(- F i 


(2.3) 


with the velocity distribution iq subject to the incompressibility constraint V • tq = 0. 
The linearised perturbation of the electromagnetic force F i = (J\Bo + JqBi) x e 2 
produced by the flow results from the interaction of the current perturbation J\ with 
the base magnetic field Bq as well as from the interaction of the base current Jq with 
the perturbation of the normal magnetic field B\. Note that the quadratic terms in the 
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perturbation amplitudes are neglected as usual in the linear stability analysis because all 
perturbations are assumed to be infinitesimal. The perturbation of the pressure gradient 
VA is eliminated from equation O by applying the curl operator. Then the flow is 
governed by the z-component of the resulting equation 


pdttoi — — (Ji • VPq + Jo ■ VPi), 


(2.4) 


where wi = e- • V x v\ is the normal component of vorticity. The associated current 
perturbation is governed by Ohm’s law for a moving medium 

Ji = a(Ei + vi x e z B 0 ). (2.5) 


A time-dependent perturbation of the magnetic field induces a rotational electric field 
according to the first Maxwell equation V x Ei = —dtB i. Note that in highly resistive 
media this induction effect is usually negligible, which corresponds to the so-called quasi- 
stationarv or inductionless approximation commonly used in the liquid-metal MHD 
( Robert j 19671 ). In order to keep our model general and applicable also to well conducting 
fluids, we forgo this approximation here. 

The charge conservation V • J\ = 0 is satisfied by introducing the electric stream 
function hi, which, as shown below, is directly related to the scalar magnetic potential 
at the current sheet. It allows us to represent the electric current in the plane of sheet as 


J i = — e z x V/ii. (2.6) 

Note that a truly 2D (z-independent) current distribution would generate a 2D magnetic 
field with z-component only which would then coincide with the current stream function 
introduced above. The current sheet considered here, however, produces a 3D magnetic 
field whose normal component at the sheet is related with the electric stream function by 
the induction-type equation which is obtained as follows. Substituting this into equation 
(E3 and taking the z-component of the curl of the resulting equation, we obtain 

d t Bi+v l -VB 0 = a~ 1 V 2 hi. (2.7) 

The magnetic field perturbation in the surrounding space is sought as 

B\ = —po V>Zq (2-8) 

using the scalar magnetic potential iq, which is governed by the Laplace equation 

V 2 tfq = 0. (2.9) 

The surface current density is related to the jump of the tangential magnetic field over the 
current sheet by Ampere’s circuital law: J i = Aj-e 2 x [_Bi] s , where = 47r x ICR' H/m 
is the vacuum permeability. Combining this current sheet condition with equations m 
and (12.81) . we obtain the boundary condition 


^U = ^iUo = ^i, (2.10) 

which relates the scalar magnetic potential with the electric stream function. It is 
important to note that the problem is defined in terms of the linear current density 
Jo which is the relevant electromagnetic parameter in the current sheet approximation. 
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3. Linear stability analysis 

Owing to the y-invariance of the base state, perturbations can be sought in the normal 
mode form 

{v 1 ,p 1 ,h 1 }{x,t) = {v,p,h}(x)e' yt+zky , 

where 7 is a generally complex growth rate and k is a real wavenumber. Then equations 
(El and El for the flow and current perturbations take the form 


1 D 2 k v = -r 1 Jok 2 (B + ^), (3.1) 

7 B = a~ 1 D 2 h + J 0 ^~, (3.2) 

2-7T X 


where v is the 2 -component of the velocity perturbation, which is further referred to as 
the transverse velocity, B is the normal ( 2 ) component of the induced magnetic field at 
the sheet, and D k = e x -^+ike y is a spectral counterpart of the nabla operator acting on 
the mode with the wavenumber k. The y-component of velocity perturbation u has been 
eliminated from (13.11) by using the incompressibility constraint D k ■ v = v' + iku = 0 , 
while the z-component of velocity is absent owing to the undeformability of the sheet. 

The boundary conditions which follow from the vanishing of the 2 -component of 
velocity and electric current at the fixed insulating edge (x = 0) and far from the edge 
(2 —> 00 ) are 

v = h = 0. (3.3) 

These boundary conditions imply ( v , h) ~ x at the edge. It means that the terms ~ 1 / 2 , 
which result from the logarithmic singularity of the base magnetic field EU> in equations 
El and El, are expected to be regular as x —>• 0 : h/x —> h'{ 0) and v/x —>■ v'{ 0). Then 
the finite longitudinal current perturbation at the edge h'( 0) = J y (0) would produce a 
perturbation of the magnetic field with a logarithmic singularity B ~ In 2 similar to 
that of the base magnetic field ( 12 . 11 ) . According to equations (13.11) and (13.21) . which 
reduce to (v, h)" ~ In x for 2 —^ 0 , logarithmic singularity produces a higher order small 
perturbation (v, h) ~ x 2 In 2 . This confirms the regularity of h 1 and {/at the edge. 

To find B required in equation El we need to solve equation El in the space 
surrounding the sheet. This can conveniently be done in the cylindrical coordinates with 
the axis aligned along the edge and the poloidal angle 9 measured from the plane of the 
sheet so that 2 = r cos0, z = rsind and 'iq (x, t) = P[r,9)e lt+lky . Then equation (12.91) 
takes the form 


rd r (rd r P) + d^P - ( rk) 2 P = 0 . 


Boundary condition (12.101) . which now reads as P 


= - W 


0=0 


0=277 


(3.4) 

= 5 / 1 ( 2 ), and the 


reflection symmetry of the problem P(r,6) = —\P(r,2ir — 9), suggest a solution in the 
form 


1 Q 1 00 

$(r, 9) = 2 H r ) cos 2 “ 2 ^ c «^«( r ) sin(ri6>), (3.5) 

~ n— 1 

where c n = sin(n0) cos | d 9 = 4ti 8 2 ' 1 _ 1 . Substituting this expression into equation (13.4L 
we obtain a sequence of ODE’s 


r{rP' n )' — ( n 2 + ( rk) 2 )P n = r{rh')' — (2 2 + ( rk) 2 )h , (3.6) 


which define the harmonics of the scalar magnetic potential x P n in terms of the electric 
stream function h. For (13.51) to be single valued at r = 0 , P n (0) = 0 is required, 
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which serves as a boundary condition for equation (13.61) . Then the normal magnetic 
field perturbation in equation (13.11) is given by 


B(x) 


d e & 


6=0 


MO 

2irx 


OO 

A nc n d' n {x). 

n =1 


(3.7) 


The problem is converted into the dimensionless form by using Vm = (fiocr) 1 and 




(3.8) 


as the characteristic velocity and time scales. The former corresponds to the magnetic 
diffusion speed, whereas the latter, when represented in terms of the characteristic 
mag netic flux density Bq = Jo, coincides with the so-called magne tic respon se time 
( Robertsl Il967 ), also known as the magnetic damping time ( Davidson 200 ill , which 
typically appears in the inductionless (Pm = 0) limit. The length scale can conveniently 
be chosen as 

t- = A. (3.9) 

where k and A are the wavenumber and wave length of perturbation, respectively. The 
electric stream function and the scalar magnetic potential are both scaled with Jo/k. 
Then equations (13.21) and (13.11) take the following dimensionless form 


(3.10) 

7 S \B = Dlh+~, (3.11) 

where the dimensionless wavenumber k = 1 and S> = tx/io Jo \J = A" is the Lundquist 
number based on the wave length A. This is the only parameter of the problem, which 
defines the Alfven speed va = Jo \j relative to that of the magnetic diffusion v m 
introduced above. A small Sa corresponds to the quasi-stationary limit in which the 
induction effect represented by the l.h.s. term of equation (13.111) becomes negligible. In the 
limit of vanishing Lundquist number (Sa = 0), the growth rate 7 becomes independent 
of S>. In the opposite limit (Sa 1), the magnetic diffusion represented by the first 
term on the RHS of equation (13.111) is expected to become negligible. In this case, 
which corresponds to an ideally conducting liquid, v can be substituted from equation 
(13.111) into equation (13.101) . Then 7 2 S| emerges as the only parameter (eigenvalue) of the 
reduced problem. Consequently, in this limit, we expect 7 ~ SJ" 1 . The respective physical 
instability development time is 


7"m 

7 



TA 

27t va ’ 


(3.12) 


which is the characteristic Alfven time for the length scale (13.91) . 


4. Numerical method 

The eigenvalue problem posed by equations (13.101) . (13.111) and (13.61) was solved nu¬ 
merically using a Chebyshev collocation method with Chebyshev-Lobatto nodes r] m = 
cos (?n. 7 r/(M + 1)), m = 0, • • • , M + 1, and the coordinate transform 77 = ^x+i 
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Figure 2. Growth rate of the leading quasi-stationary (Sx = 0) mode computed with various 
number of collocation points M and Fourier modes N versus the mapping parameter a 
(a); growth rates of the three leading modes versus the Lundquist number computed with 
M = 64, N = 4096, a = 0.3 (b). 


maps the semi-infinite domain x = r E [0, oo) onto 77 = [—1,1] using parameter 
a and transforms the differentiation operator ^ into ^a(l — ?y) 2 ^. Equations were 
approximated at the internal collocation points x m corresponding to m = 1 . ■ ■ ■ M and 
the boundary conditions were applied at xm +1 = 0 and xq = 00 (Boyd 2013 1. This 
eliminates the potential edge singularity of the normal magnetic field from the discretizied 
problem. The problem was reduced to a standard matrix eigenvalue problem as follows. 
First, the Fourier series (13.51) was truncated at the length N. Second, the normal magnetic 
field B defined by equation (13.71) was expressed in terms of the electric stream function 
h by inverting N matrices of size M x M representing the l.h.s. operators in (13.61) for 
n = 1, • • • ,N Fourier modes, and then substituted into equations (13.101) and (13.111) . 
Eventually, inverting the matrix representations of the l.h.s operators of the resulting 
equations, we obtained a standard matrix eigenvalue problem of size 2 M x 2 M for the 
growth rate 7 and the unknown vector (v, h) at the internal collocation points. In the 
quasi-stationary limit, which corresponds to Sx = 0 , the problem can be simplified further 
by expressing h from equation (13.111) in terms of v and then substituting it into equation 
(13.101) . This results in a standard matrix eigenvalue problem of size M xM for the growth 
rate 7 and the vector of v. In the ideally conducting limit, which corresponds to S\ O- 1, 
the problem is reduced to a standard matrix eigenvalue problem of the same size M x M 
for the eigenvalue y 2 S 2 and the vector h as outlined in the previous section. Matrices were 
inverted and eigenvalue problem solved using the LU and QR factorization algorithms 
from the standard linear algebra software library LAPACK. 


5. Results 

We start with the quasi-stationary limit, which corresponds to Sx = 0, and first verify 
the accuracy of our numerical method. The highest growth rate 71 computed in this limit 
with various numbers of collocation points M and Fourier modes N is plotted in figure 
[2][a) against the mapping parameter a. First, note that 71 is real and positive, which 
means a monotonic instability. The same holds also for the two subsequent eigenvalues 
shown in table [T] which are several times smaller than 71 . Second, the growth rate is 
seen to be nearly constant over a wide range of the mapping parameter a. For M = 64 
collocation points, a noticeable variation, which is a purely numerical effect, appears only 
for a < 10 ~ 2 and a > 10. The growth rate computed with M = 64 is nearly constant 
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M 

N 

7i 

72 

73 

64 

64 

0.032960 

0.014264 

0.0089024 

256 

256 

0.033375 

0.014653 

0.0092855 

1024 

1024 

0.033480 

0.014752 

0.0093836 

64 

1024 

0.033480 

0.014752 

0.0093836 

64 

2048 

0.033498 

0.014769 

0.0094000 

64 

4096 

0.033506 

0.014777 

0.0094082 


Table 1. Growth rates of the three leading quasi-stationary (S\ = 0) modes computed with 
various number collocation points M, Fourier modes N and a = 0.3. 


almost over three decades of a. As seen in figure [21 the range of constant 71 extends 
over more than five decades of a when M ^ 256. In the following, we fix a = 0.3, 
which is close to the centre of this range. Then M = 64 produces 71 coinciding up 
four decimal places with the value computed with M = 1024, which is a typically fast 
convergence of the Chebyshev collocation approximation. The convergence, however, is 
much slower with respect to the number of Fourier modes N. As seen in tableQ] N > 1024 
is required to compute the leading eigenvalue 71 « 0.0335 with 4 d.p. Convergence is 
presumably slowed down by the edge singularity, which is discussed below, as well as by 
the discontinuity of the scalar magnetic potential over the current sheet resulting from the 
boundary condition (12.101) . This discontinuity limits the Fourier series convergence rate 
to algebraic while the convergence rate o f the Chebyshev collocation approximation is 
typically exponential ( Gottlieb fc Orsza51l97?t ). Further, we use N = 4096 and M = 64. 

Growth rates of the three leading modes computed using the full dynamical model 
(Ena bud are plotted in figure [21 b) versus the Lundquist number Sa- For Sa -C 1, 
the quasi-stationary limit considered above is obviously recovered. As Sa increases, the 
development of instability is slowed down by the increasing magnetic diffusion time. For 
Sa > 10 , the growth rate is seen to reduce as 7 ~ S^ 1 , which confirms the asymptotics 
inferred in f[3J The asymptotic growth rate of the dominant 7 ~ 0.325/Sa is seen to agree 
well with the prediction of the reduced dynamical model for ideally conducting liquid 
described in the previous section. 

Amplitude distributions of the three fastest growing modes in both poorly and ideally 
conducting limits are plotted versus the distance from the edge in figure [3] The fastest 
growing mode, whose streamlines are shown in figure [4ja,b), is seen to have a single row 
of recirculation cells along the edge. Figure^ a,b) shows that each subsequent mode has 
one more row of cells, which are separated by the zero crossings of v. The perturbation 
of the normal magnetic field in the quasi-stationary limit, which is shown in figure Ele), 
has a characteristic logarithmic singularity similar to that of the base field (12.11) . This 
singularity is presumably behind the slow convergence noted above. As seen in figure 
Hf), this singularity vanishes in the ideally conducting limit. For a negligible magnetic 
diffusion, equation (13.111) at the edge (x —» 0) reduces to 7S|13(0) = v'(0)/2tt, which 
implies a finite normal magnetic field B( 0) proportional to the longitudinal velocity 
u( 0) = i'D'(O) along the edge. But this, in turn, means a singularity in the longitudinal 
current density, which is seen in figure Efd) to increase as ~ a; -1 / 2 towards the edge. 
Such a singular current distribution is typical at the ed ges o f perfectl y conducting sheets , 
where it generates a purely tangential magnetic field ( Priede et al. 20061 : Priedelfioi ll) . 
The current singularity is obviously smoothed out by finite magnetic diffusion, which 
thus produces a normal magnetic field through the sheet with a logarithmic singularity 
at the edge. The associated spatial amplitude distribution of the scalar magnetic potential 
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Figure 3. Amplitude distributions of the transverse velocity v (a,b), the longitudinal current 
density —h' (c,d) and the normal magnetic field B (e,f) for the three fastest growing modes in 
the quasi-stationary (a,c,e) and ideally conducting (b,d,f) limits. The amplitudes are normalized 
by the condition v'(0) = — 1. The normal magnetic field in the ideally conducting case is rescaled 
with (27T7 Sx) -1 so that B( 0) = {/(O). 


^{x, z) may be seen in figure [5] to have a symmetric discontinuity over the sheet with a 
continuous normal derivative which defines the normal magnetic field through the sheet. 
In the ideally conducting limit shown in figure [5jb) , the perturbation is entirely due to 
the advection of the magnetic field and, thus, more localized at the edge. In the quasi- 
stationary limit shown in figure [5j a), the perturbation is spread out over the sheet by 
the magnetic diffusion. 

As may be seen in figure [3j the increased current density at the edge is produced 
by the flow towards the edge. It means that the flow effectively compresses the current 
lines where it is directed towards the edge, while the opposite is the case where the flow 
is directed away from the edge. As this happens between the vortices, the associated 
























Figure 4. Streamlines (a,d), electric current lines (b,e) and the normal magnetic field isolines 
(c,f) of the fastest growing mode in the quasi-stationary (a,b,c) and ideally conducting (d,e,f) 
limits. 



Figure 5. Amplitude distributions of the scalar magnetic potential perturbation IF in the 
(*, 2 )-plane at the edge of sheet for the fastest growing mode in the quasi-stationary (a) and 
ideally conducting (b) limits. 


pattern of current lines in figure |UJd,e) is shifted by a quarter wave length relative to the 
streamline pattern in figure U[a,d). 

The compression of the current lines occurs as follows. According to equation (12.71) . 
a flow in a non-uniform magnetic field induces electric current in order to conserve the 
magnetic flux carried by the liquid. Since the magnetic field drops off with the distance 
from the edge, the flow towards the edge locally reduces the magnetic flux density, which 
may be seen in figure [d^f) for the perfectly conducting limit. This, in turn, implies an 
increase in the current density at the edge. The current density at the edge has to raise in 
order to compensate the magnetic field generated by the rest of the sheet. The increased 
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local current density interacting with the base magnetic field enhances the pinch force. 
This effect outweighs the advection of weaker magnetic field towards the edge. As a result, 
the initial flow perturbation towards the edge is enhanced, thus giving rise to instability 
in ideally conducting liquid. 

Additional amplification of the instability results from the magnetic diffusion, which, 
as discussed above, smooths out the current singularity, so producing a logarithmic sin¬ 
gularity in the normal magnetic field at the edge. As seen in figure[2te), this perturbation 
is positive and, thus, enhancing the base held at the edge. Therefore, the growth rate 
of instability increases as the magnetic diffusion becomes stronger and the Lundquist 
number respectively smaller. As seen in hgure[2jb), the growth rate saturates at Sa < 10, 
which means that the instability becomes dominated by the magnetic diffusion. In this 
range of Sa, the instability grows on the magnetic response time scale (13.81) . which by a 
factor of S^ 1 exceeds the respective Alfven time scale (13.121) . The latter is seen to apply 
only to Sa > 10 . 


6. Summary and conclusions 

In this study, we considered a pinch-type instability in a semi-infinite planar sheet of 
an inviscid incompressible liquid with a straight rigid edge carrying a uniform tangential 
electric current. The electromagnetic pinch force resulting from the interaction of the 
electric current with its own magnetic in this model is balanced by the pressure gradient 
in quiescent liquid. This equilibrium state was found to be inherently unstable with 
respect to the internal flow perturbation caused by a row of counter-rotating vortices 
along the edge. The vortices compress the electric current lines where the flow is directed 
towards the edge and disperse them in the opposite case. The compression and dispersion 
of the current lines occurs as a weaker magnetic field is advected by the flow towards the 
edge and a stronger held is carried away. The interaction of the current perturbation with 
the base magnetic held produces a pinch force that amplifies initial how perturbation. 
In an ideally conducting liquid, this part of the pinch force dominates over the opposing 
one which results from the interaction of the magnetic flux perturbation with the base 
current. 

The role of the latter changes in the resistive limit, where it becomes destabilizing 
and dominant. The difference in the instability mechanism is due to the current and the 
magnetic held distributions at the edge. In an ideally conducting liquid, perturbation 
of the magnetic held remains finite at the edge, whereas the associated current density 
increases as ~ x _1 / 2 becoming unbounded at the edge. In a resistive liquid, this current 
singularity is smoothed out by the magnetic diffusion, which thus produces a logarithmic 
singularity in the magnetic held. The latter enhances the base magnetic held at the edge. 
This gives rise to a destabilizing perturbation of the pinch force which dominates in 
the highly resistive liquid over the interaction of the current perturbation with the base 
magnetic held. 

Different instability mechanisms lead to disparate instability development times in the 
highly resistive and well conducting huids. In the former, the instability develops on 
the magnetic response time scale, which depends on the conductivity. In the latter, 
the instability develops on the much shorter Alfven time scale. This is because the 
perturbation of the magnetic held in poorly conducting huid depends on the product of 
conductivity and velocity, whereas in well conducting huid the magnetic held is effectively 
frozen in and, thus, its perturbation depends directly on the displacement. Therefore 
these two instability development times are generic and, thus, expected to hold also in 




12 


J. Priede 


other geometries with fixed boundaries regardless of singularities in the magnetic field or 
the electric current distributions. 

It follows from our model that the relevant electromagnetic parameter of the problem 
is the linear curren t densit y. In the eut ectic alloy of GalnSn with p = 6.4 x 10 3 kg/m 3 and 
a = 3.3 x 10 6 S/m lWeber et al. ( 201 3) carrying linear current density of Jo = lkA/cm, 
which is comparable to the respective quantity for the m = 1 instability m ode in 


the cylinder with the radiu s of 5 cm and the critical current of 2.9 kA (iRiidiger et al 


1201 lL ISeilmaver et al. 20121 . characteristic instability development time is r m / 7 i ft 
30p/(cr/iQ Jg ) ft 4 s, where 71 ft 0.0335 is the dimensionless growth rate of the most 
unstable resistive mode. Note that in a Li-Te liquid_metal battery, where the areal current 
density can reach up to jo = 7 A/cm 2 ( Kim et al. 20131) . the linear current density 
Jo = lkA/cm corresponds to a layer of thickness S ~ lm. The resp ect i ve instabi lity 
devel opment time in Li (<r = 3.3 x 10 6 S/m, p = 0.5 x 10 3 kg/m 3 (jMiiller fc Buhlen 


20011) is then by an order of magnitude shorter than that in GalnSn. Comparing the 

. . . c2 

instability development time in GalnSn with the viscous damping time t„ = due 
to the kinematic viscosity v in the respective Hele-Shaw flow (Batchel or 119671) . we can 
estimate critical Hartmann number Ha c = poJoS\/(T / pu = 12/yi ft 20. This value 

is comparable with H a„ ft 25 fo r the non-axisymmetric (to =1) instability mode in the 
cylindrical geometry (Rudiger et al. 120111 ). and it i s also not far from Ha c ft 40 for the 
axisymmetric mode in the annular geometry Priedel (2015:). For Ha > 10, a more adequate 


estimate ma y be provided by the so-called Hartmann damping time th = fpHa 


( Sommeria & Moreau 19821 ). which is based on the Hartmann rather than Poiseuille 
velocity profile, as in the Hele-Show flow, and yields Ha c = 2/71 ft 60. For liquid metals, 
which are relatively poor conductors characterized by a low magnetic Prandtl number 
Pm = p,o(Jv ~ 10 -5 — 10 -6 , all these critical Hartmann numbers ~ 10 correspond to 
the Lundquist number Sa = HaPrn 1 ^ 2 < 0.1, which means a highly resistive mode of 
instability. 

In conclusion, note that in a viscous fluid, the th resh old of instability, which is defined 
by 7 = 0, depends only the Hartmann number ( Riidiger et al. 2011). This is because 
the magnetic Prandtl number vanishes from equations (13.1011 and (13.111) together with 
the Lundquist num ber similarly to the thermal Prandtl number in the case of Rayleigh- 
Benard instability ( Chandrasekhar! 196ll) . But it is important to stress that the actual 
mechanism of the instability and, in particular, its growth rate depend on the Lundquist 
number. Also non-linear evolution of the instability depends not on only on the Hartmann 
number b ut also on the Lundquist (magnetic Prandtl) number as recently demonstrated 
bv lHerreman et al. ( 2015 ). 
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